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Abstract 

The 0(a^njTpCA,F) terms to the massive gluonic operator matrix elements are calcu- 
lated for general values of the Mellin variable N. These twist- 2 matrix elements occur as 
transition functions in the variable flavor number scheme at NNLO. The calculation uses 
sum-representations in generalized hyper geometric series turning into harmonic sums. The 
analytic continuation to complex values of N is provided. 



1 Introduction 



Heavy quark contributions to the deep inelastic scattering structure functions play a crucial role 
in the QCD analyses to determine the parton distribution functions and the strong coupling 
constant a s (M§) in a consistent manner, cf. [T]. The heavy flavor corrections were calculated 
at NLO in semianalytic form in [2"fl To avoid contributions of higher twist the analysis has to 
be restricted to large enough values of Q 2 . It has been shown in [I] that for Q 2 >10 m 2 , with 
m the heavy quark mass, the heavy flavor contributions to the structure function F 2 (x, Q 2 ) are 
rather accurately described using the asymptotic representation in which all power corrections 
oc (m 2 /Q 2 ) k ,k G N + are neglected. In this case the heavy flavor Wilson coefficients can be 
calculated analytically. They are given by convolutions of massive operator matrix elements 
(OMEs) and the massless Wilson coefficients, cf. Ref. [HE]. The massless Wilson coefficients 
are known to 3-loop order [5J. At NLO the massive OMEs were calculated in [4^ITl-fl"2"] in the 
unpolarized and polarized case, including the 0(a 2 e) contributions, and in [T3] for transversity. 
The heavy flavor corrections for charged current reactions are available at one loop and in the 
asymptotic case at two-loops [T4"],ll~5]. 

At 3-loop order a series of moments has been calculated for all massive OMEs for N = 
2. ..10(14) contributing in the fixed and variable flavor scheme, [5]. The 3-loop heavy flavor 
corrections to F L (x,Q 2 ) in the asymptotic case were calculated in [T6]. First results for general 
values of N have been obtained for the OMEs with operator insertions on the quark lines 
in case for the color factors rifT F CA,F [E] and 3-loop ladder topologies [18]. First T f Ca,f- 
contributions at general N were calculated in [19] for two heavy quark lines carrying the same 
mass. Furthermore, the moments iV = 2, 4, 6 in case of the OMEs contributing to the structure 
function F 2 (x,Q 2 ) with two different heavy quark masses were computed in [TnH^Ti] - In all the 
above cases the massive OMEs are calculated for external massless partons which are on-shell. 
The case of massive on-shell external lines has been treated in [21] recently. 

In the present paper the 3-loop corrections of 0(nfT F CA,F) to the massive OMEs with local 
operator insertions on the gluonic lines, A gq Q and A ggj Q, at general values of N are calculated. 
Together with the corresponding terms with the insertions on the quark lines, [T7J, these contri- 
butions complete all terms corresponding to the case of one massless and one massive fermion 
line at 3-loop order. These matrix elements contribute to the transition functions needed to 
describe the parton densities in the variable flavor number scheme (VFNS). In this scheme it is 
possible to define heavy quark distribution functions assuming that there exists only one heavy 
quark and all other quarks can be dealt with as massless in the sense of an effective field theory 
approach. These distributions can be used for effective calculations in some processes at hadron 
colliders. The picture holds to 2- loop orders. Starting with the 3-loop corrections, [13120], dia- 
grams containing quarks of two different masses contribute even to the universal corrections. 
Since m 2 /m 2 w 1/10 is not a small number, the original VFNS-picture does not necessarily hold 
in practice. Here we deal with the 0{njT F CA,F) contributions which are in accordance with the 
VFNS. In Section 2 the main formalism is lined out. The calculation is performed in D = 4 + e 
dimensions and uses representations in terms of generalized hypergeometric functions. They 
lead to multiple sum representations, which are solved using modern summation technologies 

1 A fast and precise numerical implementation in Mellin space has been given in [3]. 
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encoded in the package Sigma [22J. The results of the calculation are given in Section 3 both in 
Mellin-iV and in x-space, and Section 4 contains the conclusions. 



2 Parton distribution functions in the VFNS 

The neutral current Born cross section of unpolarized deep inelastic scattering (DIS) is given 
by [23] 



dxdy 



2ira 2 
xyQ 2 



M 2 ( 2 M 2 \ y 2 



T 2 {x,Q 7 



+xy(l - y)J c 3 (x,Q 2 

neglecting lepton mass contributions. Here x and y denote the Bjorken variables and — q 2 = 
Q 2 = xyS, with q 2 the 4-momentum transfer. The structure functions J~i(x, Q 2 ) contain electro- 
weak effects due to Z-boson exchange and differ for lepton and anti-lepton-nucleon scattering, 
cf. [23], and 



i . 



(2) 



In the limit M| 3> Q 2 the electromagnetic terms in (CQ) dominate and only the two structure 
functions Fi 2 (x,Q 2 ) contribute, with 



2xF 1 (x, Q 2 ) = F 2 (x, Q 2 ) - F L (x, Q 2 



(3) 



where F^ is the longitudinal structure function. Both structure functions contain light and 
heavy quark contributions. The y-dependence of the differential scattering cross section is used 
to separate the structure functions [24j and allows precise measurements of the structure function 
F 2 (x, Q 2 ). In the twist- 2 approximation, referring to the fixed flavor number scheme, they are 
given by 



F 2 {x, Q 2 ,n f ) = F 2 m=0 (x, Q 2 , n f ) + F~ (x, Q 2 , n f , m) 



(4) 



Here F™ °(x, Q 2 ) denotes the well-known massless contribution and the massive contribution in 
the presence of a single massive quark reads [5] 



2,0 



\x,Q 2 ,n f ,m) = ^2eU L^ q fn f , <8> f k (x,fi 2 ,n f ) + f-i:(x,fi 2 ,n f ) 

k=i i fn [i j 

~ P c f Q 2 m 2 \ 2 ~<- / Q 2 m 2 \ 2 . 

+ L 2, g [rif, ^2,^2)® S(a;, /*,"/) + L 2g \n f , ^ , J ® G(x, [i , n f ) 



+4 



p S / Q 2 m 2 
2,5 I n /> 9 » ~ 2 



) S(x, Li 2 ,rif) 

/ Q2 m 2\ 

+ H l 9 in/, ^ , J G(x, /i 2 , n/) 



(5) 
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Here fk{x, /i 2 , nf), /i 2 , nf), G(x, /i 2 , rif) denote the kth quark, singlet-quark, and gluon den- 
sities, respectively with 



n f , /i 2 ) = ^[/*(a:, n f , /x 2 ) + fo(x, n f , fi 2 )} . 



(6) 



k=l 



The Wilson coefficients (n / , Q 2 / m? , m 2 / /i 2 ) and Lf (n/, Q 2 /m 2 , m 2 / n 2 ) have been calcu- 
lated completely for general values of iV in [TT] . 

The renormalization group implies the following representation for the set of (rif + 1) (mass- 
less) parton densities expressed in terms of n/ parton densities j8] : 

/kfa + 1, /x 2 , m 2 , iV) + fr(n f + 1, /i 2 , m 2 , N) = < S Q f n,, ^, N) ■ [f k (n f , /i 2 , iV) 

+/> / ,/i 2 ,A0] 



rn 



/ g (n / + l,/i 2 ,m 2 ,iV) + /Q(n / + l,/i 2 ,m 2 ,iV) = Ag! ( n f , ^, iV ) • E(n,, /i 2 , iV) 



2 



-KV Q n/ ,-^,JV •£(/*/, /i 2 , AO, (7) 



+A Qg n/.-^.JV \-G{n f ^\N) . (8) 



Here /q(/q) are the heavy quark densities. The flavor singlet, non-singlet and gluon densities 
for (rif + 1) flavors are given by 



E(n f + l,ti 2 ,m*,N) 



2 



N 



E(n /; /i 2 ,iV) 



+ 



// 



PV> -o . ^ + A Qg n f 



m 



' 2 ' 



N 



G(n f ,fi 2 ,N) 



A(n f + l,ti 2 ,m 2 ,N) = f k (n f + 1, fi 2 , N) + ftfn f + 1, /i 2 , m 2 , AT) 

1 -£(n/ + l,/x 2 ,m 2 ,#) 



n/ + 



(9) 
(10) 



G(n f + l,fM\m\N) = A gq , Q [n f , ^, N ■ S(n/, /x 2 , AT) 



Any relation between the (nf + 1)- and n^-parton density can only contain universal, i.e. process- 
independent, quantities. 

Note that the new parton densities depend on the renormalized heavy quark mass m 2 . As 
outlined above, the corresponding relations for the operator matrix elements depend on the 
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mass-renormalization scheme, with m = m(a s (/j, 2 )) in the MS scheme, which we will apply 
below. These equations describe the transition of one heavy quark becoming light at the time 
referring to the scale fi 2 . 

The matching scales /i 2 are often chosen as /i 2 = m 2 . The comparison of the results in 
complete calculations to those in which flavor thresholds are matched in the VFNS allows in 
principle to determine the relevant matching scale. In an analysis of the various deep-inelastic 
structure function sum rules [25] it has been shown that the scale fi 2 turns out to be significantly 
different of m 2 . This is not unexpected since mass effects do not turn into the behaviour of the 
massless case close to the production threshold. 

The resummation of large logs, as being performed in the VFNS, has to be performed at very 
high scales. As has been shown in [26] this is not the case in the kinematic range at HERA. A 
smooth transition from the threshold region to asymptotic scales has been proposed in terms of 
the BMSN-scheme gj, 

Ff(x,Q 2 ,n f = 4) = F?' FFNS (x,Q 2 ,n f = 3) + Ff ,asymp (x,Q 2 ,n f = 4) - Ff> asymp (x, Q 2 , n f = 3) ,(12) 



which is found to be in excellent agreement with the HERA data [27]. There is a series of 
other proposals to match between the threshold and asymtotic region [28H30] . partly with a 
faster transition to the massless case. Here precise data on -F 2 cc (x, Q 2 ) are helpful to distinguish 
between different descriptions. We would like to mention that a correct treatment of the heavy 
flavor corrections is of instrumental importance in the QCD analysis of the complete structure 
functions F 2 (x,Q 2 ), which has been measured to a precision of 0(1%) [31J. 



3 The 0{oL Z s n f Tl) contribut ions to Agg^Q and A gq ,Q 



The OMEs A gq ^Q and A gg ^Q are expectation values (j\O g \j), i,j — q,g of the gluonic operator 



O 



g,m,...,n N 



2i N - 2 SSp[F flia D fl2 ...D m _ 1 F« N ] - trace terms . (13) 

between massless on-shell external states. The corresponding massive OMEs A gqt Q, A g9t Q were 
calculated to 0(a 2 ) in [8] and including also terms linear in e in [9] correcting the previous result. 

The renormalized expressions A gq ^Q and A gg; Q to 0(al) were derived in [5]. In the MS scheme 
with the heavy quark mass m on-shell they are given by : 



A 



(3), MS 

gq,Q 



(o) 



24 



~(o)a(o) _J_ Ly(0) 
igq '19 1 V >qq 



7 <°> + 10 A, + 24/3 ,q ) A),Q }ln*{ — )+^ 67^0, 



m 



m \ '99 



7 <f - 4/3 - 6/? 0lQ ) 



+ 7™(7r S + 7«' PS -7fi> + 2 A , Q ) l^(-) 



I 1 



(2),PS 

l Qq 



(2) 
l 99,Q 



Iqq 



^(iff-iS , + <A + «flw)+iff(- 



7<»» + 12/3„, Q + 10ft 



ft,, 



/m 



(0)^(2) 

L gg,Q 



_(2),PS _(2),NS 



a. 



Qq 



qq,Q 



(0)/- 

Tot s>3 (_ (0)~(0) 



24 



+ 

'99 '9S 



7^-7^ + 10/30 



' .'-/// 



9<?,Q 



(14) 



5 



.4 



(3), MS 



99,' 



Here Smf^ are expansion coefficients of the unrenormalized mass, /^g are coefficients of the 
/3-functions (including mass effects), (k is the Riemann-^ function with k G N\{0, 1}, a\j ,a[j are 
two loop contributions to order e° and e l respectively, and 7^, %j are the anomalous dimensions, 
and quantities with a hat or a tilde are defined by 



48 



+4/3 2 + 14/3 , Q /3 + 12/3 2 q )/3 ,q| ln 3 (^) + ± J'yW ( T « + (1 - n,)^)) 
+7i?7^ + H'J^Q ~ VfStfo + 2A),q] + 4[A + /3 1iQ ]/3 ,q 



<S? '99 ' '99 ^ u >v '99 

m 2 \ 1 



ln^ — +- 



/i- 



16 



5 '99 ° n f a gq,Qiqg 



16a® g (2A) + 3A),g) 



+87S )o S + 8 ^°Xq + 7S ) 7£ ) C 2 (7^ - 7$ + 6/5 + 4n / /3 , Q + 6/3 , Q 
+4/30,0(2(7^ + 2/?o) (2^0 + 3/3 ,q) I In(^) + 2(2/3 + 3/^)^ 



+ 



(0) * (0) > 

ri /7?9 a gg,Q '99 tt Q9 '99 + 4g \'w '99 z \- zn f + L \P°> C 

A),qC3 



-6/3 J + 

*(0)> 

7VC2 

16 



12 



[Po,Q ~ 2A>W + 2[/?o + 6/3 ,q]A),q - 4# 
7«+ 7^+^(7^ -27^-2^-2^) + 



(-1) 



(2) 
Z 99,Q 



+245mS 0) /3 , Q + SSm^g + C 2 /3o,q/3o + + Sm? [f3 , Q 6m 

+8m? (^) 7 W) + 2A W 7£ ) + 4/3 ,q/3o + 8/3 2 Q ) - 2<W, 0) /? ,q + a^g . (15) 



,(0) 
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/ = /(n / + l)-/(n / ), 



(16) 



see Ref. [5]. The unreormalized OME A^g also receives contributions from the vacuum polar- 
ization insertions on the external lines 



flj(p 2 , m 2 , /x 2 , a 2 ) = [-^ 2 + p,p v ) & (k \ 

n« = n«(o, m 2 ,/i 2 ) 



2-2 2\ 

p ,m ,fi) 



k-1 



such that 



2(3) _ J(3),1PI _ A(3) _ i(2),lPIA(l) _ „ 5(1) fr(2) , J(l) fr(l)fr(l) 



,(3,0) 
l 99," 



,(3,1) 



,(3,2) 



+ 



a gg,Q + a gg,Q 



+ a 



(3) 

gg,Q ' 



All contributions to (j!4)15p but the constant terms a^'g are known [il[7H9|[T2|[32] . In particular, 
all the logarithmic contributions have already been obtained for general values of the Mellin 
variable N, 



(17) 
(18) 

(19) 
(20) 
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In the following we calculate the contributions O^a^nfTpCp^) to the massive gluonic OMEs. 
The Feynman diagrams are generated by QGRAF [33] and the extension allowing to include local 
operators [5]. The color-algebra is performed using |35j. For a large part of the calculation 
we use FORM [36]. The momentum integrals are performed introducing a Feynman parame- 
terization. The Feynman parameter integrals are then rewritten in terms of hypergeometric 
functions (2-^1,3-^2), which are represented in terms of absolutely convergent series. The re- 
sulting sums, which may still contain finite sums due to binomial expansions, are then pro- 
cessed applying the symbolic summation technology, which is encoded in the package Sigma [22] 
and making use of a large number of algorithms for processing multi sums using the package 
EvaluateMultiSums [37j[38]. Additionally it is very useful to reduce such sums to a smaller 
number of 'key sums', by synchronization of the summation ranges and algebraic reduction of 
the summands. This step helped to reduce the size of the terms from 2GByte to 7.6MByte 
and the number of sums from 2419 to 29. The algorithms for this step are implemented in the 
package SumProduction [39]. Details of the corresponding technique are described in [37] [38]. 
The corresponding expressions have simplified using mutual relations and methods applicable 
to the respective classes of sums encoded in the package HarmonicSums [3D]. The results for the 
individual diagrams have been checked comparing to the moments obtained in [5] using the code 
MATAD [U]. The constant contributions afjQ, j — q,g to (I14|ll5p read : 

(3WT| 2 f 16(iV 2 + iV + 2) (1 2 o 



16 (8iV 3 + 13iV 2 + 27iV + 16) 

27(iV - 1)N(N + l) 2 
32 (35iV 4 + 97iV 3 + 178iV 2 + 180iV + 70) 
27(iV - 1)N(N + l) 3 

32 (1138iV 5 + 4237iV 4 + 8861iV 3 + 11668iV 2 + 8236iV + 2276) 



(3C 2 + S 2 + S 2 ) 
Si 



2A3(N - 1)N(N + i; 



1 



(21) 



(3),n/T| _ 2 J r 1 

a gg,Q ~ n f I F\ C A 



(N-l)(N + 2) 



4P X 8P 2 

+ 0/ 



27iV 2 (iV + l) 2 1 729iV 3 (iV + l) 



160 448 Pi 

+-^-(N -1)(N + 2)( 2 Si - -?-(N -1){N + 2)C 3 5! + 



27 v yv ^ 27 v yv ^ 729iV 4 (iV + l) 4 

2P 4 j. 56 (3iV 4 + 6iV 3 + 13iV 2 + ION + 16) 4P 5 

rC2 H , -t \ C3 — n „, r0/ , r , -,x 9 ^2 



27iV 2 (iV + l) 2 ^ 27iV(iV + l) ^ 27N 2 (N+l)< 



(N-l)(N + 2) 



2 



U2(N 2 + N + 2Y 3 _ 16P 6 
27iV 2 (iV + l) 2 1 27iV 3 (iV + l) 3 1 



32P 7 i 6 (iV2 + iV + 2) 2 i 6 (iV2 + iV + 2) 2 

H O A 7"9 / A 7" 1 1\9 S2&1 H O A 7"9 / A 7" I 1 N 2 ^2^1 



U7V 4 (iV + l) 4 3iV 2 (iV + l) 2 ^ 3iV 2 (iV+l) 

32P 8 16P 9 448 (N 2 + N + 2) 2 16P 10 



>C2 H nAT9/Ar , 1N9 C3 + nsr „„ , 7\3$ 



243iV 5 (iV + 1) 5 9iV 3 (iV + l) 3 ^ 9iV 2 (iV + l) 2 *° 9iV 3 (iV + r 

160 (N 2 + N + 2f 



2 



27iV 2 (iV + l) 2 



-S, 



(22) 
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where the polynomials Pj are given by 



Pi 

P-2 



Pi 
P 5 
P G 

Py 



P9 
Pio 



16N 5 + 41N 4 + 2N 3 + 47N 2 + 70N + 32 (23) 
6944iV 8 + 26480iV 7 + 23321iV 6 - 15103iV 5 - 39319iV 4 - 27001iV 3 - 11178iV 2 
-2016iV + 864 (24) 
4809iV 10 + 24045iV 9 - 182720iV 8 - 854414iV 7 - 1522031iV 6 - 1472927iV 5 
-758234iV 4 - 126080iV 3 - 1152iV 2 - 50688iV - 24192 (25) 
3iV 6 + 9iV 5 + 307iV 4 + 599iV 3 + 746iV 2 + 448iV + 96 (26) 
40iV 6 + 112iV 5 - 3iV 4 - 166iV 3 - 301iV 2 - 210iV - 96 (27) 
44iV 6 + 123iV 5 + 386iV 4 + 543iV 3 + 520iV 2 + 248iV + 24 (28) 
205iV 8 + 856iV 7 + 3169iV 6 + 6484iV 5 + 731(W 4 + 4722iV 3 + 1534iV 2 
+48iV - 72 (29) 
1976iV 10 + 9385iV 9 + 24088iV 8 + 38989iV 7 + 50214iV 6 + 53872iV 5 + 35219iV 4 
+6890iV 3 - 4233iV 2 - 2844iV - 756 
14iV 6 + 33iV 5 + 59iV 4 + 39iV 3 + 55iV 2 + 20iV - 12 



AN b + 3N b - 50iV 4 - 129iV 3 - lOOiV 2 - 56iV - 24 . 



(30) 
(31) 
(32) 



Here S^s == Sb,a(N) = J2 n =i s ig n (b) n 5'a(-/V)/'^' b '; S$ = 1 denote the harmonic sums [12] which 
only occur as single harmonic sums in the present calculation. 
It is convenient to express the renormalized OMEs A g j t Q, j 

quark mass in the MS scheme, cf. 0. The OMEs A®£ ,7 % and A?}% ,7 % read : 



q, g also referring to the heavy 



99,< 



A 



(3), MS 

gq,Q,C F T^n f 



J3),n f Tj,,MS 
^99,Q 



C F n f T\ 



, 2) 32(N 2 + N + 2) 
F * 9(N -1)N(N + 1) 



+ 



16 (iV 2 + iV + 2) ( a . 32 (8iV 3 + 13iV 2 + 27iV + 16) 

[b l + D 2 J H ^7— , ; no ^1 



3(JV-l)iV(JV+l) 



9(N -1)N(N + 1) 



+ 



32 (19iV 4 + 81iV 3 + 86iV 2 + 80iV + 38) 
27(N -1)N(N +1) 3 

32 (N 2 + N + 2) 



(S* 3 + 3S2S1 + 2^3 

.Si 



27(N -1)N(N + 

32 (8iV 3 + 13iV 2 + 27iV + 16) 

27 {N - l)iV(iV + l) 2 
64 (4iV 4 + 4iV 3 + 23iV 2 + 25iV + 8 
27(iV - 1)N(N + l) 3 

64 (197iV 5 + 824iV 4 + 1540iV 3 + 1961iV 2 
243(iV - 1)N(N + l) 4 

64 (N 2 + N + 2) 2 
9(N - 1)N 2 (N + 1) 2 (N + 2) 



128(N 2 + N+ 1) 64 
27(iV - 1)N{N + 1){N + 2) ~ 27 1 



24C 3 ) 



1388iV + 394) 




+C A 



In 3 



m 



(33) 



8 



16 . o f m 2 



-C^ln 2 ™ + [C A 



\N - l)(7V + 2) 



16Pi2 
817V 2 (7V+1) 2 1 



1 



4P 



13 



"97V 4 (7V+1) 4 37V 3 (7V+1) 3 



(JV — 1)(7V + 2) 
32Pi4 



4Pn 
817V 3 (7V+ l) 3 

16 (TV 2 + 7V + 2^ 
7V 2 (7V + 1) 2 



Si 2 - ^ : 



(7V-l)(7V + 2) 




8P, 



^1 2 

+^(7V-1)(7V + 2)C 3 S 1 



27 



4P 



277V 2 (7V + 1) 2 2 
1 



^ ID q 

7297V 3 (7V +1) 3 1 

1024 (TV 2 + A/" + 1) 
7297V 4 (7V+ l) 4 277V(7V + 1) ^ 3 



277V 2 (7V + l) 2 
2Pi 7 



+C 



'(7V-l)(7V + 2) 
32P 19 



64 (TV 2 + TV + 2f 



9N 2 (N + 1) 2 
64P 



1 4 N 

-g^i 3 - 8C 3 + 



32P 



21 



~27TV 3 (TV+ l) 3 1 81TV 4 (TV + 1) 4 1 243TV 5 (TV + l) 5 
32P 22 



3TV 3 (TV + l) 3 
with the polynomials 



5 2 



Pn = 

P 2 = 

Pl3 = 

Pu = 

Pl5 = 

Pie = 

Pn = 

Pl8 = 

Pl9 = 

P20 — 

P21 = 



297iV 8 + 1188iV 7 + 640iV 6 - 2094iV 5 - 1193TV 4 + 2874TV 3 + 5008iV 2 
+3360iV + 864 

136iV 6 + 390iV 5 + 19TV 4 - 552iV 3 - 947iV 2 - 630iV - 288 

15TV 10 + 75TV 9 - 48TV 8 - 866TV 7 - 2985iV 6 - 6305iV 5 - 8206iV 4 - 7656TV 

-4648TV 2 - 1600iV - 288 

r5 



(34) 



(35) 



(36) 



(37) 
(38) 



(39) 
(40) 
(41) 



5TV 5 + 52iV 4 + 109TV 3 + 90TV 2 + 48TV + 16 
4TV 5 + 17TV 4 + 14TV 3 + 71iV 2 + 70iV + 32 

3008iV 8 + 11600iV 7 + 9197iV 6 - 10255iV 5 - 27739iV 4 - 24745TV 3 - 124747V 2 
-20167V + 864 (42) 
41857V 10 + 209257V 9 + 18927V 8 - 1171187V 7 - 2221517V 6 - 1768637V 5 - 414467V 4 
+223047V 3 - 12967V 2 - 184327V - 6912 (43) 
167V 6 + 527V 5 - 37V 4 - 1067V 3 - 2777V 2 - 2107V - 96 (44) 
107V 6 + 307V 5 + 1097V 4 + 1687V 3 + 1557V 2 + 767V + 12 (45) 
387V 8 + 2067V 7 + 9627V 6 + 22467V 5 + 25097V 4 + 15427V 3 + 5097V 2 + 247V - 36 (46) 
1237V 12 + 7387V 11 + 6917V 10 - 35267V 9 - 145217V 8 - 294587V 7 - 391897V 6 
-376727V 5 - 219207V 4 - 39147V 3 + 28567V 2 + 18727V + 432 



P 22 = 27V 6 + 47V 5 + TV 4 — 107V 3 - 57V 2 - 47V - 4 . 



(47) 
(48) 
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As has been noted before [5J, the above results are free of (2, which is common to all massive 
OMEs, and hence is a particular feature of representing also the mass in the MS scheme. Fur- 
thermore we note, that the In 2 (m 2 //i 2 )-contribution to ^ggQc F T 2 n is particularly simple, while 

the corresponding contribution to _ 2 vanishes. 

As a by-product of the calculation we obtain the corresponding contributions to the anoma- 
lous dimensions from the single pole term 1/e resp. the linear logarithmic contribution, cf. ( fl4lfT5l) . 



2 ( U(N 2 + N + 2) , a c ^ 128(8iV 3 + 13iV 2 + 27iV + 16) c 
UfTFCF [~3(N-1)N(N+1) [Sl + S " } + 9(N-1)N(N+1) 2 Sl 



128 (AN 4 + AN 3 + 23N 2 + 25N + 8) 
9(JV- 1)N(N+ l) 3 



(49) 



Igg — n f l F Lsj 



+n f T 2 F C F 



32P23 c 8P24 



27(N -1)N 2 (N +1) 2 (N + 2) 27(N - l)iV 3 (iV + l) 3 (iV + 2) 

64 (i y 2 + N + 2 ) 2 _ 
3(iV - 1)N 2 (N + 1) 2 (N + 2) 1 1 2; 

lfiP„, 

(50) 



128P 25 16P 26 



9(iV- l)iV 3 (iV + l) 3 (iV + 2) 27(iV- l)iV 4 (iV + l) 4 (iV + 2) 



where 



P 23 = 8iV 6 + 24iV 5 - 19iV 4 - 78iV 3 - 253iV 2 - 210iV - 96 (51) 
P 24 = 87iV 8 + 348iV 7 + 848iV 6 + 1326iV 5 + 2609iV 4 + 3414iV 3 + 2632iV 2 + 1088iV 

+ 192 (52) 

p 25 = AN 6 + 3iV 5 - 50iV 4 - 129iV 3 - lOOiV 2 - 56iV - 24 , (53) 
P 26 = 33A^ 10 + 165A^ 9 + 256A^ 8 - 542A^ 7 - 3287A^ 6 - 8783A^ 5 - 11074A^ 4 - 9624A^ 3 

-5960iV 2 - 2112iV - 288 . (54) 

Eqs. (I49II50I) confirm previous results in [32] by a first direct diagrammatic calculation, here in 
the massive case. 

The leading singlet eigenvalue for the gluonic anomalous dimensions 7^ - , j = q,g in form of 

(2),«? (0) 

iT ) + ^ (55) 

199 n f 

has been calculated in |13] for the leading rif contribution, oc n 2 . We also confirm this result by 
a direct massive calculation. 

Usually the calculation in iV-space is being performed multiplying the massive OMEs and 
the parton distributions analytically^], cf. e.g. [44] . The corresponding analytic continuations of 
harmonic sums up to weight w=8 are given in [35] . Only a single numerical contour integral 
around the singularities has to be performed, allowing for very fast implementations. 



2 For Mellin-space representations of a wide class of parton densities see [13] . 
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The OMEs 



can also be given in x-space directly for codes operating in x-space only. 



They are given by : 



A (3),n f 1f,MS M 



gg,Q ^ > 



C F n f Tp< 



32a; 64 64 



9 



9x 9 



m 



16.r 32 32'. , 
1 I Hi 

3 3a; 3 



+ 

+ 
+ 



/256a; 320 320 
+ 



V 9 

32a; 



27 



+ 



9;r 
64 



9 



27a; 27 



6 4m 



608s 2176 2176 

Hi H 1 

27 27a; 27 

256a; 320 320 



/256x 128 128 
\~27~ ~ 27x + ~W ' ' 

24064 24064 1 



27 27.r 

256a; 



27 



H? 



9 



512 512\ 12608a; 
~9x + ~9~)^ + 243 



243a; 



243 



7 



(56) 



n f T j 



64a; 2 64a; 



64 



27 



27 



64a; 128 „ XTT 64 256 

1 (1 + x)H H 1 

9 9 V ; 9 27a; 



27(z- 1)+ 
In 3 



128 64 

~ri + 27x 



— o 

m 



256a;' 
27 



16 

C F 5(1 - x) In 2 

3 



+ 



32 44 
+ f(l + ^-p(l 



C A 
16 



608a; 2 16 , 

144C 2 -85 )x 

27 81 V S ; 



(144C 2 + 149) 



832x 2 16a; 

+ 17 



27 



500 



27 
256 



I He 



832x 2 208a; 



(1 + z)H 0)1 - 



27 9 
2176 



+ 



9 v " ' '~ u,i 81(x- 1) + 
256a; 2 688a; 592 



9 



9 



9 



176 832 
IT + 27x 

224 |+c F 

64a; 2 



27a; 



Hi 
32 



(l + x)H 3 



-f (12C 2 + 5)a; 



64 



(12C 2 



+ 



-41) )HcH 

64a; 2 



512a; 2 128a; 128 512\ 
9 3" + ~3~ + ~9x) 



256 

)HiH + — (l+x)H ,iH 



16a; + 16 + 



64 



64 

3x 



H? 



90 64 64 

y 5(l-x) + -x 2 (18C 2 -7) + -x(3( 2 



+3C 3 -28) + -(6C 2 + 3C 3 + 10) + 



+ 



128x 2 64a; 256 512 
9 + ~3 3~~~9x 



64a;' 



256 



416a; 736 896 . T1 

1 Hi 

3 3 9a; ' 



H ,i - — (1 + x)H 0A i + 64(1 + x)H ,i,i 



3904 
^7x 



In 



m 



/512a; 2 512a; 



512 



V 27 



27 27(x - 1). 



+ 



1024 512' 



27 



27x 



11 



.'2048a; 2 512a; 1024,., XTT 512 204 
+C F ( 3 ( nfy + — — (1 + a;)H ~ 



+ 



27 9 9 

208a; 2 812a; 320 



9 



27x 



125 
"81 



(l + x)Hg 



81 



+ 



64 



27(x - 1) 
128 



81 
H + 



+ 



81 



H 



8624a; 2 8 , _ . 16 , , 

^ir-8T(«C 2 -199),--(8C 2 + 19) 



416a; 2 56a; 88 416 \ TT TT /64 > 310\ 
-"sT + lT "27 + 8li HlH ° + 27 C2 -^ )5(1 - X) 



w. TT /208a; 2 20a; 44 208 \ rT , 
— (1 + aOH 0il Ho + ^— + — - — ) H 2 



27 



27 81a; 



416 



x 2 (9C 2 + H3) - — (2088(2 - 864C 3 + 6055) - — a;(2601C 2 - 864C 3 - 4883) 



729 



+ 



729" ' 729' 

8624a; 2 2600a; 872 4592 \ /832a; 2 2144a; 2120 416 ' 



243 



81 



+ 



125 
~1?7 

+C F 



(l + x)(H 0) o,i + H ,i,i) 



81 243a; , 

24064 



V 81 
32320 



81 



81 



81a; 



+ 



729 (a; - 1)+ 729a; 



128a; 2 256a; 64 



+ 



^ H 



-|(18C 2 -1))h 2 +(- 

/128a; 2 32a; _ 32 _ 128 \ 
' V 27 + ~9~ ~ ~9 ~~ 27x) 

64 \ 
(225C 2 - IO8C3 - 1292) H 



81 81 SI 

128a; 2 32a; 32 12 



9 



+ Y + 27a^J H 1 H + -(l + a;)Ho )1 H 



H?H + ( -||(18C2 - l)x 2 - ||(333C 2 - IO8C3 - 410)x 



243 



+ 



4352a; 2 320a; 704 896 ., . , . 

1 1 HxHq 

81 9 9 81a; ' 



/512a; 2 2560a; 1408 256 \ 128,, ri , T1 ... 

+ I — — + — — + ~ — I H ,iH - — (1 + x) (H ,o,iH + Ho,i,iHo) 



+ 
+ 



27 27 

/256a; 2 64 



27 



27a; 



64. t , 256 \ tt3 / 



V 81 

/1024 >2 1312\ P/i 



9 

1472a; 2 64. 1472,.,, 

(x — 1) H Hf 

81 9 V J 81x 1 



V 9 



-cl 



64 64 
—a; (63C| + 145C 2 - 120C 3 + 1720) + ^(^(2 



64 



-IO8C3 - H65) - ^ ( 63C| - 215C 2 - 30 + C 3 - 1675 



^(18C 2 -D, 2 



+ ^[(3C 2 + 44)a;-(3C 2 + 80)] 



128(18C 2 -163) \ 



243a; 



/1408a; 2 128 r , rt> or7X 
+ ( -gj— + ^[(9(2 + 37)a; + (9C 2 - 71)] 



396 
ttx~ 



Ho,i + 



512a; 2 
27 



2560a; 
27 



1408 256 

+ — IH 



27 



/256a; 2 1664a; 1664 256 
+ 



27.i Z" ' ' 1 + V 27 27 27 27. v 



)h 0i i,i + ^(1 + x) ^Ho.o.0,1 



+Ho,o,i,i _ 2110,1,1,1 + 



79744 
729a; 



(57) 
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with the harmonic polylogarithms = Hg(x) over the alphabet 21 = {0,1,-1} [15]. They 
can be expressed in terms of elementary functions and the Nielsen integrals [47J : H (x) = 
ln(x),Hi(x) = -ln(l -ar),H ,i(s) = Li 2 (x), H 0A i(x) = Li 3 (x), H ,i,i(x) = Si i2 (x), H ,o,o,i(^) = 
Li 4 (x),H 0A i i i(x) = S 2j2 (z) and H ,i,i,i(a;) = Si, 3 (ar), with 

S n , P {x) = \ > / -!Lhi<^%)h?Q.-xy), (58) 

(n -ly.pl Jo y 

hi n (x) = S n - 1>n (x) . (59) 

Here Li n (x) denotes the polylogarithm. All higher functions but ^^(x) can be reduced to poly- 
logarithms by the argument relation x — > (1 — x). Numerical implementations of the functions 
S n>p (x) were given in [18] . 

At small values of x the functions A^^ f Q F ' MS (x) are singular as oc 1/x, or in A^-space like 
oc 1/(N — 1), unlike the quarkonic contributions given in [TTj with a leading pole oc 1/N. One 
notices that the number of functions needed in x-space to express A^£?q f ' MS is larger than in 
iV-space, as has been found also in other analyses, cf. [UHUHS], requesting very careful numeric 
implement at ions . 



4 Conclusions 

We have calculated the contributions 0(alrifTpCA,F) to the massive OMEs with local operator 
insertions on gluonic lines and veritces at general values of the Mellin variable N. These matrix 
elements are needed to describe the transition functions in the VFNS. In the calculation repre- 
sentations of the Feynman diagrams by generalized hypergeometric functions play an essential 
role. They allow the e-expansion into nested sums, which can be solved using modern summation 
technologies. The number of these sums is very large, although their structures exhibit simi- 
larities. One may synchronize these sums, leading to a low number, however, with voluminous 
intermediate terms. The solution of the latter sums turns out to be more economic. The final 
results in N space can be expressed by rational functions in N and single harmonic sums up to 
Ss(N). We also derived the corresponding x-space results, which have a more involved structure 
and depend on six Nielsen integrals. 
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